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Abstract. We have obtained medium resolution spectra 
of the neutron star X-ray transient Aql X— 1 during qui- 
escence. We determine the spectral type of the compan- 
ion star to be Kl and also estimate its rotation speed 
to be 62i|gkms _1 . By measuring the width of the Ha 
emission line profile of the accretion disc, we estimate the 
binary inclination to ~ 50° and also estimate the semi- 
amplitude of the companion star's radial velocity curve to 
be ~l70kms _1 . 

Key words: binaries: general - stars: fundamental pa- 
rameters, stars: individual: Aql X-l - stars: neutron 



1. Introduction 

Aquila X-1(=V1333 Aquilae) is a soft X-ray transient 
source that shows X-ray bursts (Koyama et al. 1981; Cz- 
erny, Czerny & Grindlay 1987), thereby indicating that 
the compact object is a neutron star. From quiescent ob- 
servationss the companion star has been identified to be a 
V=19.2 K-type star. The optical counterpart brightens by 
~2-5 magnitudes during X-ray outbursts which is inter- 
preted as reprocessing of radiation in the accretion disk 
(Thorstensen, Charles & Bowyer 1978, hereafter TCB; 
Canizares, McClintock & Grindlay 1979; Charles et al. 
1980; van Paradijs et al. 1980). 

Attempts to find the orbital period have revealed many 
variations but no firm period. Watson (1976) reported an 
unconfirmed 1.3 day X-ray periodicity during the 1975 
outburst. Chevalier & Ilovaisky (1991) have obtained an 
18.97 hr periodicity from optical photometry during its 
active state, which they interpret as being the orbital pe- 
riod. 

TCB observed the optical counterpart in its low state. 
They found the only significant feature in the spectrum 
was the Mg6 5 175 A blend in absorption; no Balmer lines 
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were evident. They concluded that the spectrum is char- 
acteristic of a KO-3 V star. Shahbaz et al. (1996) also 
obtained low resolution spectra of Aql X-l in quiescence. 
They too found the Mgb absorption blend but also Ha 
and the Paschen series in emission. They found by fitting 
the continuum spectrum that it could be best described 
by a K5 star. 

Aql X-l is known to undergo regular X-ray and op- 
tical ourbursts on a timescale of ~1 year (Kaluzienski et 
al. 1977; Priedhorsky & Terrell 1984; Charles et al. 1980) 
much more frequently than the other neutron star tran- 
sient Cen X-4 (McClintock & Remillard 1990). Recently 
the RXTE All Sky Monitor showed Aql X-l to be have 
undergone an X-ray outburst between late January and 
early March 1997 (Levine & Thomas 1997). Ilovaisky & 
Chevalier (1997) reported that Aql X-l was optically in 
quiescence by 30 March 1997. In this letter we report on 
our medium resolution spectra of Aql X-l obtained in 
May 1997, when the source was in quiescence. We were 
able to determine the spectral type of the companion star 
and also its rotation speed. 

2. Observations and data reduction 

We obtained intermediate resolution spectra of Aql X-l on 
13 May 1997 using the 3.9-m Anglo-Australian Telescope 
at Siding Spring, Australia. A Tek 1024 2 CCD attached 
to the RGO spectrograph was used during all the obser- 
vations. The 600V grating was centered at 5190A giving 
a dispersion of 1.55 A pixel -1 and a spectral resolution of 
3.6A (FWHM=165 km s -1 at Ha). We took 10 spectra 
of Aql X-l each with an exposure time of 1800 sees start- 
ing at 14:58 UT with interspaced Cu-Ar arc spectra for 
wavelength calibration. Template field stars of a variety 
of spectral types were also observed. 

The data reduction and analysis was performed us- 
ing the Starlink FIGARO package, the PAMELA routines of 
K. Home and the molly package of T. R. Marsh. Removal 
of the individual bias signal was achieved through sub- 
traction of the mean overscan level on each frame. This 
was acceptable since an examination of the bias frames 
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showed no significant structure. Small scale pixel-to-pixel 
sensitivity variations were removed by multiplying by a 
flat-field frame prepared from observations of a tungsten 
lamp. One-dimensional spectra were extracted using the 
optimal algorithm of Home (1986), and calibration of the 
wavelength scale was achieved using 4th order polynomial 
fits giving an rms scatter of l/12th of a pixel (0.13A). 



Table 1. Optimal Subtraction of the Companion Star 



Star 


Sp. Type 


X 2 (DOF=907) 


vsin i 


/ 


HR 805 


G8m 


794.2 


64 


1.05±0.03 


HR 794 


KOm 


779.5 


62 


1.09±0.04 


HR 807 


Kim 


776.6 


62 


0.94±0.03 


HR 822 


K2m 


927.2 


66 


0.77±0.03 


HR 872 


K3m 


900.6 


66 


0.60±0.02 


HR 871 


K4m 


1040.4 


36 


0.55±0.02 


HR 851 


K5m 


861.6 


50 


0.70±0.02 



3. The spectral type and rotational broadening of 
the companion star 

We determine the spectral type of the companion star 
by minimizing the residuals after subtracting different 
template star spectra from the Doppler-corrected aver- 
age spectrum. This method is sensitive to the rotational 
broadening v sin i and the fractional contribution of the 
companion star to the total flux (/; 1-f is the "veiling 
factor"). 

First we determined the velocity shift of the individual 
spectra of Aql X-l with respect to each template star spec- 
trum by the method of cross-correlation (Tonry & Davis 
1979). The Aql X-l spectra were then interpolated onto a 
logarithmic wavelength scale (pixel size 80 kms -1 ) using 
a sin x/x interpolation scheme to minimize data smooth- 
ing (Stover et al. 1980). No significant velocity variations 
were observed (see section 4), so we averaged the Aql X-l 
spectra to the rest frame of the template star. 

In order to determine the rotational broadening v sin i 
we follow the procedure described by Marsh, Robinson & 
Wood (1994). Basically we subtracted a constant repre- 
senting the fraction of light from the template star, multi- 
plied by a rotationally broadened version of that template 
star. We broadened the template star spectrum from to 
200 km s _1 in steps of 1 km s _1 using the Gray rotation 
profile (Gray 1976). We then performed an optimal sub- 
traction between the broadened template and averaged 
Aql X-l spectra. The optimal subtraction routine adjusts 
the constant to minimize the residual scatter between the 
spectra. The scatter is measured by carrying out the sub- 
traction and then computing the x 2 between this and a 
smoothed version of itself. The constant, /, represents the 
fraction of light arising from the template spectrum, i.e. 
the secondary star. The optimal values of vsin i and / 



are obtained by minimising x 2 (see Table 1). The above 
analysis was performed in the spectral range 5120-6710 
A excluding the Ha emission line and the interstellar fea- 
ture at 5190 A. A linear limb-darkening coefficient of 0.76 
was used (Al-Naimiy 1978) appropriate for 6000 A and 
an effective temperature of 4500 K (typical for a K star). 
We also performed the above analysis using zero and full 
limb-darkening and only obtained a change in v sin i of at 
most 4 kms -1 . 

From Table 1 it can be seen that the minimum x 2 oc- 
curs at spectral type Kl with a vsin i of 62±|g kms _1 (l- 
a) and the companion star contributing about 94% to the 
observed flux at ~ 6000 A. Fig. 1 shows the results of the 
optimal subtraction and the x 2 ~~ v sin i plot is shown in 
Fig. 2. All the template stars used in the above analysis 
were of luminosity class III. The effects of using differ- 
ent luminosity class spectra in the analysis was estimated 
by comparing the absorption line flux density of stars of 
different spectral type and class. We find using template 
stars of spectral type Kl, K4 and K7 and luminosity class 
v and III, the flux density of the stars of different luminos- 
ity class (but the same spectral type) are only different at 
the 1 percent level. Therefore, the value we obtain for / 
in the above analysis is not dependant on the luminosity 
class of the template star used. 

In order to verify the validity of the minimum in the 
X 2 - v sin i plots we performed the same analysis on a 
template star. The broadening in the template star spec- 
trum is dominated by the instrumental broadening. To 
this template star we added noise to produce a spectrum 
of comparable quality to our Aql X-l data. The above 
broadening and optimal subtraction procedure was then 
repeated, The results are shown in Fig. 2, where the min- 
imum x 2 occurs when no broadening is used. This is ex- 
actly as expected since the template star spectrum must 
be dominated by instrumental broadening. Therefore we 
are confident that the minimum x 2 for the Aql X-l spec- 
trum is real, although the errors in the v sin i determina- 
tion are large. 

4. The Aql X I spectrum 

In Fig. 1 we show the variance-weighted average of the 
Aql X-l spectra, which has a signal-to-noise ratio of 
about 25 in the continuum. The most noticeable fea- 
ture is double peaked Ha emission which has a broad 
base (FWZI^2000kms _1 ). We tried to see if there was 
any velocity variation in the absorption lines by cross- 
correlating the individual Aql X-l spectra with a tem- 
plate star (HR 807) , but found no significant variation over 
the 5 hr baseline. Therefore we quote the mean heliocen- 
tric systemic velocity obtained from the averaged Aql X- 
1 spectrum of 57±21kms~ 1 (the template star used, HR 
807, had a heliocentric velocity of 15±16kms _1 ). Again 
from the averaged Aql X-l spectrum the equivalent width 
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Fig. 1. The results of the optimal subtraction. From top to bottom: the variance-weighted average spectrum of Aql X-l, the 
template Kim star (HR 907) broadened by 62kms~ 1 , the residual spectrum of Aql X-l after subtracting the template star 
times /=0.94. The spectra have been normalized and shifted vertically for clarity. The right inset shows a close-up of the 
smoothed double-peaked Ha emission line profile in the Aql- X-l spectrum. The left inset shows a close-up of the Ha residual 
spectrum with FWHM=816±40 kms" 1 . 



of Ha is 5.3±0.3A with a mean heliocentric velocity of 
99±30kms _1 , obtained using a single Gaussian fit. 

5. Discussion 

5.1. The binary inclination 

After subtracting the spectrum of the companion star 
from the Aql X-l spectrum, the most noticable feature 
is the single peaked Ha emission line (presumably aris- 
ing from the accretion disc) with a FWHM velocity of 
816±40kms~ 1 (see Fig. 1). This is in contrast to the ob- 
served Aql X-l spectrum which shows a double-peaked 
profile, implying that the core of the double-peaked profile 
is due to the Ha absorption line arising from the secondary 
star. 

We can compare the accretion disc Ha emission line 
in Aql X-l with that in Ccn X— 4 since both have neu- 
tron star primaries and probably similar orbital periods. 
The width of the disc Ha emission line in Cen X-4 is 
687±20kms- 1 (FWHM; Casares et al. 1997) whereas it is 
broader in Aql X-l, which suggests that the binary incli- 
nation of Aql X-l may be higher than that of Cen X-4. If 
we assume Keplerian motion in the accretion disc and that 
the mass of the compact object, and orbital period of the 
two systems are comparable, then the velocities in the disc 
simply scale with sin i (Frank, King & Raine 1992). Using 
the velocity widths given above for Aql X-l and Cen X- 
4 and the binary inclination for Cen X-4 (40°; Shahbaz, 
Naylor & Charles 1993), we estimate i ~ 50° for Aql X-l. 



It should however, be noted that the assumption of Ke- 
plerian velocities at the outer edge of the disc may not hold 
(Orosz et al. 1994). For a system at such an inclination 
one might have expected to sec a double-peaked emission 
line profile, arising from an accretion disk viewed at high 
inclination. However, there are probably other sources of 
Ha emission from the system which contaminate the disc 
profile, such as from the heated face of the secondary star 
or from the bright spot. Also note the poor resolution and 
signal-to-noisc of our data, therefore we cannot rule out 
a high binary inclination on the lack of a double-peaked 
profile alone. 

5.2. The K 2 velocity 

If the orbital period is 18.97 hr and the secondary star 
fills its Roche lobe, then the density of the secondary star 
is p=0.3 gem -3 . A Kl main sequence star would have 
p=1.9 g cm -3 ., implying that the radius of the secondary 
star must be about twice that of a main sequence star 
in order for it to fill its Roche lobe. It is therefore likely 
to be evolved, similar to the secondary star in Cen X-4 
(Shahbaz, Naylor & Charles 1993). If we believe the esti- 
mate to the binary inclination to be ~ 50°, and that both 
systems have similar mass compact objects, then one can 
estimate K 2 for Aql X-l [K 2 scales with sin i c.f. the mass 
function equation). Using if 2 =146 kms~ 1 (McClintock & 
Remillard 1990) and i = 40° for Ccn X-4 and i = 50° for 
Aql X-l we estimate K 2 for Aql X-l to be ~170kms _1 . 
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Fig. 2. The X 2 - v sin i analysis plot for Aql X-l (shown as the 
histogram) using the template star HR 907. The minimum x 2 
(Xmin) is found with a rotational broadening of 62 kms -1 . The 
68% confidence level (Xmin + 1) is a l so marked as the dashed 
horizontal line. Also shown is the result for simulated data 
using a template star shown as the solid line. Note that Xmin 
occurs at zero broadening, which implies that the minimum in 
the Aql X-l x 2 ~ wsin i plot is real. The right abscissa refers 
to the template star x 2 values. 
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6. Conclusions 

Using medium resolution optical spectra, we have deter- 
mined the spectral type of the companion star in the neu- 
tron star X-ray transient system Aql X-l, to be a Kl star. 
By optimally subtracting different broadened versions of 
the companion star spectrum from the average Aql X-l 
spectrum we determine the rotational broadening of the 
companion star to be 62^20 kms -1 , with a contamination 
of ~ 6% to the observed flux at 6000A. We also estimate 
the binary inclination to be ~ 50° by measuring the width 
of the accretion disc Ha emission line, and K 2 to be ~ 170 
kms -1 . 
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